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The mass spectrum of orbitally and radially excited states of B and D mesons 
pg ■ is calculated in the framework of the relativistic quark model. The expansion 

in inverse powers of the heavy quark mass is carried out up to the first order, 

C — ' 

qv . while the light quark is treated without expansion. We find that the relativis- 

tic treatment of the light quark plays an important role. Different patterns of 
Q_i| P level inversion are discussed. The obtained masses of orbitally and radially 

O ' excited states are in accord with available experimental data and heavy quark 

symmetry relations. 
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I. INTRODUCTION 

The investigation of the properties of excited heavy-light mesons represents an interesting 
and important problem. The experimental data on the orbitally [0] and radially Q excited 
B and P meson states are becoming available now. The P and c — t factories will provide 
more accurate and comprehensive data for these states. On the other hand, the presence 
of the heavy quark in these systems considerably simplifies the theoretical description of 
heavy- light mesons. The heavy quark symmetry || , arising in the limit of infinitely heavy 
quark mass tuq, imposes strict constraints on the mass spectrum of heavy-light mesons [|J. 
In this limit the heavy quark mass and spin decouples and all the meson properties are 
determined by light degrees of freedom alone. As a result the heavy quark spin degeneracy 
of the levels emerges. The light quark's spin couples with its orbital momentum, resulting 
for 5-wave mesons in two degenerate j = 1/2 states, corresponding to 1~ and 0~ []. For 
P-wave mesons two degenerate j = 3/2 states, (1 + , 2 + ), and two degenerate j = 1/2 states, 
(0 + , 1 + ), arise. The heavy quark symmetry constrains also the strong decays of these states. 
The decay rates of the P-wave states in degenerate doublets should be the same. The states 
with j = 1/2 are expected to be wide, because they decay in S'-wave, while j = 3/2 mesons 
are narrow, because they decay in P-wave. Since the real b and c quarks are not infinitely 
heavy, the corrections in inverse powers of the heavy quark mass turn out to be important. 
These corrections break down the degeneracy of the levels. The heavy quark effective theory 
(HQET) (see e.g. || and references therein) provides a systematic method for treating I/tuq 
corrections. However, in order to obtain quantitative predictions it is necessary to combine 
it with some dynamical nonperturbative approaches. 

Many different approaches have been used for the calculation of orbital and radial exci- 
tations of heavy-light mesons @-§ . However, almost in all of them the expansion in inverse 
powers not only of the heavy quark mass (ttiq) but also in inverse powers of the light quark 
mass (m q ) is carried out. The estimates of the light quark velocity in these mesons show that 
the light quark is highly relativistic (v/c ~ 0.7 -r 0.8). Thus the nonrelativistic approxima- 
tion is not adequate for the light quark and one cannot guarantee the numerical accuracy of 
the expansion in inverse powers of the light quark mass. In this paper we use the relativistic 
quark model ||T0"|-|12]| for the calculation of the masses of orbitally and radially excited B and 



D mesons without employing the expansion in l/m q . Thus the light quark is treated fully 
relativistically. Concerning the heavy quark we apply the expansion in 1/rriQ up to the first 
order. Our relativistic quark model is based on the quasipotential approach in quantum 
field theory. It has been used for the calculation of heavy quarkonia mass spectrum |10| and 



electroweak decays of heavy mesons [pLiynfl . Recent nonperturbative investigations indicate 
that the confining potential cannot be simply scalar. This is in agreement with our model 
assumptions on the Lorentz structure of quark confinement. The comparison of our model 
results with model independent constraints of HQET is given in [T3I . 



The paper is organized as follows. In Sec. |I| we describe our relativistic quark model 
giving special emphasize on the role of the Lorentz-structure of quark confinement. In Sec. [Tl I 



1 We use the J notation, where J = j ± 1/2 and P are the total angular momentum and parity 
of the meson. 



we construct the quasipotential of the interaction of a light quark with a heavy antiquark. 
We use the heavy quark I/tuq expansion to simplify the construction. The light quark is 
treated relativistically First we consider the inn — > oo limit and then the corrections of 
the first order in I/ttiq. We present the predictions of our model for orbitally and radially 
excited states of D, D s , B and B s mesons. In Sec. [V] we compare our results with the heavy 
quark symmetry and other quark model predictions as well as available experimental data. 
Section [V] contains our conclusions. 

II. RELATIVISTIC QUARK MODEL 

In the quasipotential approach a meson is described by the wave function of the bound 
quark- antiquark state, which satisfies the quasipotential equation |TJ| of the Schrodinger 
type [|15[ 



b 2 (M) p 2 \ r d 3 q 



\ 2/j R 2/jLrJ r J (27r) 

where the relativistic reduced mass is 



*m(p) = / 7^4^(p, q; M)* M (q), (1) 



_ E q E Q _ M 4 - (m 2 q - m 2 Q ) 2 
^ R ~E q + E Q ~ 4M3 ' [Z) 

and E q , Eq are given by 

M 2 — m 2 n + m 2 M 2 — m 2 + m 2 n 

E, = 2 £ E Q = 2 S (3) 

q 2M ' Q 2M ' K J 

here M = E q + Eq is the meson mass, m qy Q are the masses of light and heavy quarks, and 
p is their relative momentum. In the centre of mass system the relative momentum squared 
on mass shell reads 

h 2 (M) _ [M 2 ~ K + m Q ) 2 ][M 2 - (m q - m Q ) 2 ] 
[ } AM 2 ' [ } 

The kernel V(p, q; M) in Eq. (|I]) is the quasipotential operator of the quark- antiquark 
interaction. It is constructed with the help of the off-mass-shell scattering amplitude, pro- 
jected onto the positive energy states. An important role in this construction is played 
by the Lorentz-structure of the confining quark- antiquark interaction in the meson. In the 
literature there is no consent on this item. For a long time the scalar confining kernel has 



been considered to be the most appropriate one [16]. The main argument in favour of this 
choice is based on the nature of the heavy quark spin-orbit potential. The scalar potential 
gives a vanishing long-range magnetic contribution, which is in agreement with the picture 
of quark confinement of Ref. [|17 |, and allows to get the fine structure for heavy quarkonia 



in accord with experimental data. However, the calculations of electroweak decay rates 
of heavy mesons with a scalar confining potential alone yield results which are in worse 
agreement with data than for a vector potential [|T^,|T^]. The radiative Ml-transitions in 
quarkonia such as e. g. J/tp —>■ r\,rj are the most sensitive for the Lorentz-structure of the 
confining potential. The relativistic corrections for these decays arising from vector and 



scalar potentials have different signs fl8|]19|| . In particular, as it has been shown in Ref. |19 
agreement with experiments for these decays can be achieved only for a mixture of vector 
and scalar potentials. In this context, it is worth remarking, that the recent study of the 
qq interaction in the Wilson loop approach pOf indicates that it cannot be considered as 



simply a scalar. Moreover, the found structure of spin-independent relativistic corrections is 



not compatible with a scalar potential. A similar conclusion has been obtained in Ref. |21 
on the basis of a Foldy-Wouthuysen reduction of the full Coulomb gauge Hamiltonian of 
QCD. There, the Lorentz-structure of confinement has been found to be of vector nature. 
The scalar character of spin splittings in heavy quarkonia in this approach is dynamically 
generated through the interaction with collective gluonic degrees of freedom. 

All these new results are in agreement with the assumptions of the relativistic quark 
model approach [|I9| , |T0|| . Indeed, in constructing the quasipotential of quark- ant iquark in- 



teraction we have assumed that the effective interaction is the sum of the usual one-gluon 
exchange term with the mixture of long-range vector and scalar linear confining potentials, 
where the vector confining potential contains the Pauli interaction. The quasipotential is 
then defined by 



+V c v onf (k)^ q T Q ,, + kL (k) Jt^uoHf), ( 5 ) 

where as is the QCD coupling constant, D^ is the gluon propagator in the Coulomb gauge 
£°°(k) = -g, D*(k) = -^ (& - *?f\ , D« = i7« = 0, (6) 

and k = p — q; 7^ and u(p) are the Dirac matrices and spinors 

u\p) 



\ 2e(p) ' CTP ]X ' 



crp 



with e(p) = y/p 2 + m 2 . The effective long-range vector vertex is given by 

r M (k) = 7M + ^^. ( 8 ) 

where k is the Pauli interaction constant characterizing the anomalous chromomagnetic 
moment of quarks. Vector and scalar confining potentials in the nonrelativistic limit reduce 
to 

C(r) = (l-e)(A- + B), 

V c s oni (r)=e(Ar + B), (9) 

reproducing 

Konf (r) = Klf (r) + Konf (r) = Ar + B, (10) 



where e is the mixing coefficient. 

The expression for the quasipotential for the heavy quarkonia, expanded in v 2 /c 2 , can 
be found in Ref. |I0| . The structure of the spin-dependent interaction is in an agreement 



with the parameterization of Eichten and Feinberg |[22| . All the parameters of our model 



like quark masses, parameters of linear confining potential A and B, mixing coefficient e and 
anomalous chromomagnetic quark moment k were fixed from the analysis of heavy quarkonia 
masses [IU| and radiative decays |TjJ. The quark masses m& = 4.88 GeV, m c = 1.55 GeV, 



m s = 0.50 GeV, m u ^ = 0.33 GeV and the parameters of the linear potential A = 0.18 
GeV 2 and B = —0.30 GeV have standard values of quark models. The value of the mixing 
coefficient of vector and scalar confining potentials e — — 1 has been determined from the 
consideration of the heavy quark expansion |13| and meson radiative decays [JHJ. Finally, 



the universal Pauli interaction constant k = —1 has been fixed from the analysis of the fine 
splitting of heavy quarkonia 3 Pj- states JTOj] . Note that the long-range magnetic contribution 
to the potential in our model is proportional to (1 + k) and thus vanishes for the chosen 
value of k = —1. 



III. QUASIPOTENTIAL OF THE INTERACTION OF A LIGHT QUARK WITH 

A HEAVY ANTIQUARK 

The expression for the quasipotential (H) can, in principle, be used for arbitrary quark 
masses. The substitution of the Dirac spinors ([?[) into (|5]) results in an extremely nonlocal 
potential in the configuration space. Clearly, it is very hard to deal with such potentials 
without any simplifying expansion. Fortunately, in the case of heavy-light mesons, one can 
carry out (following HQET) the expansion in inverse powers of the heavy quark mass tuq. 
The leading terms then follow in the limit tuq — > oo. 

A. Infinitely heavy quark limit 

In the limit tuq — > oo the heavy quark Dirac spinor has only an upper component. Thus 
only the zeroth component of Dirac matrices in ([5]) contributes, and we get 

4 4vr 

7 g ° + £-7 g °(7fc) +V^(k)\u q (q). (11) 

The resulting quasipotential is still nonlocal in configuration space. However, taking into 
account that B and D mesons are weakly bound 0, we can replace e q (p) — > E q in the Dirac 
spinors (|7|). Such simplifying substitution is widely used in quantum electrodynamics [p3} -p5 




and introduces only minor corrections of order of the ratio of the binding energy (V) to E q . 
This substitution makes the Fourier transformation of the potential ( |TT| ) local and works well 



2 The sum of constituent quark masses tuq + m q is very close to the ground state meson mass M. 



for the confining part of the potential. However, it leads to a fictitious singularity 1/r 3 at 
the origin arising from the one-gluon exchange part, which is absent in the initial potential. 
Note that this singularity is not important if the perturbation theory in l/m q is used. As we 
are not going to expand our potential in l/m q , additional analysis is necessary. The explicit 
one-gluon contribution to (O) is given by 



Uqip) - 



4 
•— < 

3 



47T 



u q {q) 




l(p) + 



m. 



k/(p) + 



2e q (p) \ 2e q (q) 

pq + i<T g [pq] 



M) + 



m n 




(12) 



The most interesting is the last term in this expression, which comes from the lower com- 
ponents of the spinors u q . In configuration space, it has an 1/r asymptotics for r — > and 
decreases as 1/r 3 for r — > oo. If we simply replace e q — > E q , then we get the asymptotics 
1/r 3 , which is correct at infinity but too singular at the origin. To cure this asymptotics, let 
us notice that if binding effects are taken into account, it is necessary to replace e q — > E q — V, 
where V is the quark interaction potential. If we consider larger distances of order of the 
hadron radius, the contributions of V are of order of the (small) binding energy and can 
be neglected. But for r — > the Coulomb singularity in V becomes important. Thus, we 
keep this contribution (up to the first order in the binding energy) in the last term in (|T2"|). 
As a result we find that the Coulomb-like one-gluon potential Vc ou i( r ) = — | v" m * n ^ s term 
should be replaced by 



Vboui(r-) = Vcoui(r) 



1 



(n-!t;)(n-! 



E q +m„ r 



(13) 



leading to the correct asymptotics both at the origin and infinity. 

The resulting local quark- ant iquark potential for rriQ — ► oo can be presented in configu- 
ration space in the following form 



VmQ^ooyT) 



E q + m q 
2E q 

+V c V on{ (r 



Vcom(r) + V cou{ (r) + 

" KonfWlP 



1 



(E q + m q y 



2m n 



+ ~ Kcoul(0 



v:L(r) - v'7 



conf( r ) 



conf 



P[Vcoul(^) 

(r)[l-(l + K)] 
-2(1 + ^^ 



LS r 



(14) 



Here the prime denotes differentiation with respect to r, L is the orbital momentum, and S q 
is the spin operator of the light quark. Note that the last term in (0) is of the same order 
as the first two terms and thus cannot be treated perturbatively. 

In the infinitely heavy quark limit the quasipotential equation ([I]) in configuration space 
becomes 



m' 



P z 



2E q 2E Q/ 



^Af(r) = V n , 



(r)*Af(r) 



(15) 



and the mass of the meson is given by M = itlq + E q 



Solving flT5|) numerically |26| we get the eigenvalues E q and the wave functions ^m- The 
obtained results are presented in Table d We see that the heavy quark spin decouples in 
the limit iriq — > oo, and thus we get the number of degenerated states in accord with the 
heavy quark symmetry prediction. 



B. 1/rriQ corrections 

The heavy quark symmetry degeneracy of states is broken by 1/rriQ corrections. The 



corrections of order I/tuq to the potential (14) arise from the lower components of heavy 
quark Dirac spinors in (|5|). The matrix elements of the spatial components of Dirac matrices 
are not zero now. Thus the contributions at first order in l/rriq come from the one-gluon- 
exchange potential and the vector confining potential, while the scalar potential gives no 
contribution at first order. The resulting 1/rriQ correction to the heavy quark potential (|14"D 
is given by the following expression 



SVi 



l/m Q \ 




IIV i 
conf ' 



>]) 



+ 3 



-S g S Q + — (S g r)(S Q r) 



AlWr) + (1 + n) 2 AV 



conf V 



SqSq 



(16) 



where S = S q + Sq is the total spin. The first three terms in (|16"D represent spin-independent 
corrections, the fourth term is responsible for the spin-orbit interaction, the fifth one is the 
tensor interaction and the last one is the spin-spin interaction. It is necessary to note that 
the confining vector interaction gives a contribution to the spin-dependent part which is 
proportional to (1 + k) or (1 + k) 2 . Thus it vanishes for the chosen value of k — — 1, while 
the confining vector contribution to the spin-independent part is non zero. 

The quasipotential at 1/rriQ order is given by the sum of V m .^^oo(r) from (fH]) and 



SVi/ mQ (r) from ([16]). By substituting it in the quasipotential equation ([[]) and treating the 
1/rriQ correction term SVi/ mQ (r) using perturbation theory, we are now able to calculate 
the mass spectrum of D, D s , B and B s mesons with the account of I/tuq corrections. The 
results of our calculations are presented in Tables 0-0. 



IV. RESULTS AND DISCUSSION 



Let us compare the obtained results with model independent predictions of HQET. In 
HQET the heavy-light meson mass M of ground-state pseudoscalar and vector mesons is 
given by 



M = \x Q + A 



Ai - 2 [J( J + 1) - 

7 



A, 



o{i/fi 



(17) 



where the parameter A represents contributions arising from terms in the HQET Lagrangian 
which are independent of the heavy quark mass. The terms Ai and A2 parametrize contri- 
butions from the kinetic energy and chromo-magnetic interaction. 

Note that the HQET heavy quark mass /xq is in principle different from our mass tuq 
by the simple reason, that HQET uses the pole heavy quark mass, while quark models use 
the constituent quark mass. Obviously, in the heavy quark symmetry limit this difference 
between constituent and pole heavy quark masses, however, disappears. In particular, we 
find the relation E q — A + 0(l/rriQ). Thus, the values of E q listed in Table | are just the 
values of the HQET parameters A. 

The structure of l/mg corrections in our model is consistent with (|iT|). However, we find 
that the parameter corresponding to Ai in our model contains not only the kinetic energy 
part but also other terms originating from spin- independent corrections in Eq. (|T6| ) 3 . A 
similar result has been found in the relativistic quark model in Ref. ||27|| . The value of this 
parameter turns out to be very sensitive to the value of the heavy quark mass. Indeed, 
assuming that the heavy quark pole mass \iq and constituent quark mass rriq are related by 

m 

mQ = HQ-^— + 0(l/tJL 2 Q ), (18) 

we see that only Ai is influenced, and it is connected to the quark model value A^ by 

Ai = A? M + \ Q . (19) 

In our model, for constituent quark masses mj, = 4.88 GeV and m c = 1.55 GeV, we find 
Ai ~ 0.85 GeV 2 for B mesons and Ai ~ 0.26 GeV 2 for D mesons. If we require the 
parameter Ai to be equal to the mean value of —(p 2 ) (in our model (p 2 ) ~ 0.25 GeV 2 for 



B and D mesons), then from ([ll|) and (|l^) we get the heavy quark masses /i;, = 4.75 GeV 
and /i c = 1.40 GeV. These values agree with values of the heavy quark pole masses used in 
HQET. 

The values of the parameter A2, which determines the hyperfine splitting, coincide in 
HQET and quark models. We find A 2 ~ 0.112 GeV 2 for B mesons and A 2 « 0.125 GeV 2 for 
D mesons. 

Heavy quark symmetry provides relations between excited states of B and D mesons, 
such as 

M Bl - M Dl = M Bal - M Dsl =M B -M D = M Bs - M Ds = m b - m c , (20) 

where M Bl = (3M Bl + 5M B2 )/8, M B = (M B + 3M#*)/4 are appropriate spin averaged P- 
and S"-wave states. We get from Tables [TT|-[V] the following values of mass splittings 

M Bl - M Dl = 3.29 GeV 



3 Note that, if we had used also an expansion in inverse powers of the light quark mass, then, in 
the static limit, we would get Ai = —{p 2 }- 

8 



fB Bl ~ M Dsl 


= 3.30 GeV 


M B -M D 


= 3.34 GeV 


M Bs - M Da 


= 3.33 GeV 



(21) 

in agreement with (p0[). There arise also the following relations between hyperfme splittings 
of levels 

777 777 

AM B = M B2 - M Bl = —AM D = -£ (M Da - M Dl ) , (22) 

m b m b 

and the same for B s and D s mesons as well as for P\ — Pq states. Our model predictions for 
these splittings are displayed in Table [VT]. 

In Tables 0"0 we compare our relativistic quark model results for heavy-light meson 
masses with the predictions of other quark models of Godfrey and Isgur |5|], Isgur ||, Eichten, 
Hill and Quigg |7| and experimental data p8l,|i"|,|2|j . All these quark models use the expansion 
in inverse powers both of the heavy rriQ and light m q quark masses for the Qq interaction 
potential. In Ref. |J some relativization of the potential has been put in by hand, such as 
relativistic smearing of coordinates and replacing the factors l/m q by l/e q (p). However, the 
resulting potential in this approach accounts only for some of the relativistic effects, while 
the others, which are of the same order of magnitude, are missing. The considerations of 
Refs. ||0] are closely related. The heavy quark expansion is extended to light (u, d, s) quarks 
and the experimental data on P wave masses of K mesons are used to obtain predictions 
for B and D mesons. 

In the paper || it is argued that the heavy quark spin P-wave multiplets with j = 1/2 
(0 + , 1 + ) and j = 3/2 (l + ,2 + ) in B and D mesons are inverted [2?J. The 2 + and 1 + states 



lie about 150 MeV below the 1 + and + states. In the limit mg — > oo, we find the same 
inversion of these multiplets in our model, but the gap between j = 1/2 and j = 3/2 states 
is smaller (~ 90 MeV for B and D mesons and ~ 70 MeV for B s and D s mesons), and 1/tuq 
corrections reduce this gap further. However, the hyperfme splittings among the states in 
these multiplets turn out to be larger than in 0. As a result, the states from the multiplets 
for D, D s and B s mesons overlap in our model, however the heavy quark spin averaged 
centres are still inverted (see Figs. 1-4). We obtain the following ordering of P states (with 
masses increasing from left to right): B meson -- 1Pi(|), IP2, lP(b LPl(§); D s meson - 
1P , 1Pi(|), 1P 2 , 1Pi(|); D and B s mesons -- 1Pi(|), 1P , 1P 2 , lPi(f). Thus only for B 
meson we get the purely inverted pattern. Note that the model predicts the ordinary 
ordering of levels. 

The results of our model agree well with available experimental data. The experimental 
values in Tables |H|-[V| for ground state and P-wave masses are taken from Refs. [^1,0. For 
the radially excited states we use the preliminary data from DELPHI 0. 



V. CONCLUSIONS 

In this paper we have presented the calculation of the mass spectra of orbitally and 
radially excited states of heavy-light mesons in the framework of the relativistic quark model. 
The main advantage of the proposed approach consists in the relativistic treatment of the 
light quarks (u,d,s). We apply only the expansion in inverse powers of the heavy quark 



(b, c) mass, which considerably simplifies calculations. The infinitely heavy quark limit as 
well as the first order l/rriq corrections are considered. Our model respects the constraints 
imposed by heavy quark symmetry on the number of levels and different splittings. 

We find that the heavy quark spin multiplets with j = 1/2 (0 + , 1 + ) and j = 3/2 (1 + , 2 + ) 
are inverted in the ttlq — > oo limit. This inversion is caused by the following reason: The con- 
fining potential contribution to the spin-orbit term in ([14]) exceeds the one-gluon exchange 
contribution. Thus the sign before the spin-orbit term is negative, and the level inversion 
emerges. However, 1/itlq corrections, which produce the hyperfine splittings of these mul- 
tiplets, are substantial. As a result the purely inverted pattern of P levels occurs only for 
the B meson. For D and B s mesons the levels from these multiplets begin to overlap. This 
effect is more pronounced in the D s meson, where the ordinary ordering is restoring. Thus 
we see that the D s meson occupies an intermediate position between heavy-light mesons and 
heavy- heavy mesons (quarkonia, B c meson). It is necessary to note, that both the presence 
of the relativistic light quark and the light-to-heavy quark mass ratio play an important role 
in the formation of level ordering patterns. The light quark determines the meson radius, 
while the m q /rriQ ratio indicates the validity of the application of the heavy quark symmetry 
limit. 

The found mass values of orbitally and radially excited heavy-light mesons are in good 
agreement with available experimental data. At present only narrow P-wave 2 + , 1 + (j = 3/2) 
levels of D, D s and B mesons have been measured. It will be very interesting to observe also 
+ , 1 + (j = 1/2) levels, which is more complicated, because these states are expected to be 
broad. This will allow to determine the ordering of P levels and to test quark dynamics in 
a heavy-light meson. We plan to apply the found meson wave functions for the calculation 
of semileptonic and nonleptonic decays of B mesons into P-wave D and D s mesons, which 
are important for the experimental observation of these states. 
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TABLES 

TABLE I. The values of E q and meson masses M in the limit wiq — ► oo (in GeV). We use the 
notation [n 3 L) for meson states, where n is the radial quantum number, and j is the total angular 
momentum of the light quark. 







D 






D s 






B 






B s 




State 


E u ,d 




M 


E s 




M 


E u ,d 




M 


E s 




M 


lis 


0.497 




2.047 


0.607 




2.157 


0.514 




5.394 


0.628 




5.508 


lip 


0.800 




2.350 


0.913 




2.463 


0.800 




5.680 


0.927 




5.807 


lip 


0.886 




2.436 


0.985 




2.535 


0.898 




5.778 


0.998 




5.878 


2iS 


0.940 




2.490 


1.047 




2.597 


0.955 




5.835 


1.063 




5.943 



TABLE II. Mass spectrum of D mesons with the account of 1/rnn corrections in comparison 
with other quark model predictions and experimental data. All masses are given in GeV. We use 
the notation (nLj) for meson states, where J is the total angular momentum of the meson. 



State 



our 



I 



experiment 



■ 



ISo 


1.875 


1.88 




lSi 


2.009 


2.04 




m 


2.459 


2.50 


2.460 


m 


2.414 


2.47 


2.415 


m 


2.501 


2.46 


2.585 


m 


2.438 


2.40 


2.565 


25 


2.579 


2.58 




25i 


2.629 


2.64 





1.8645(5) 
2.0067(5) 
2.4589(20) 
2.4222(18) 



2.637(9) 



TABLE III. Mass spectrum of D s mesons with the account of I/tuq corrections in comparison 
with other quark model predictions and experimental data. All masses are given in GeV. The same 
notation as in Table |0| is used for meson states. 



State 



our 



§ 



experiment 



■ 



ISo 


1.981 


1.98 




lSi 


2.111 


2.13 




1P 2 


2.560 


2.59 


2.561 


lPi 


2.515 


2.56 


2.526 


IPi 


2.569 


2.55 




1-Fb 


2.508 


2.48 




2S 


2.670 


2.67 




25i 


2.716 


2.73 





1.9685(6) 

2.1124(7) 

2.5735(17) 

2.535(4) 
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TABLE IV. Mass spectrum of B mesons with the account of 1/rriQ corrections in comparison 
with other quark model predictions and experimental data. All masses are given in GeV. The same 
notation as in Table |H| is used for meson states. 



State 



our 



1 



experiment 



lSb 


5.285 


5.31 






5.2792(18) 


ISi 


5.324 


5.37 






5.3248(18) 


1P 2 


5.733 


5.80 


5.715 


5.771 


5.730(9) 


m 


5.719 


5.78 


5.700 


5.759 




ipi 


5.757 


5.78 


5.875 






IPo 


5.738 


5.76 


5.870 






2S 


5.883 


5.90 








25i 


5.898 


5.93 






5.90(2) ? 



TABLE V. Mass spectrum of B s mesons with the account of 1/mQ corrections in comparison 
with other quark model predictions and experimental data. All masses are given in GeV. The same 
notation as in Table |H| is used for meson states. 



State 



our 



I 



1 



experiment 



131 



ISo 


5.375 


5.39 




ISi 


5.412 


5.45 




1P 2 


5.844 


5.88 


5.861 


IPi 


5.831 


5.86 


5.849 


IPl 


5.859 


5.86 




IP) 


5.841 


5.83 




25 


5.971 


6.27 




25i 


5.984 


6.34 





5.3693(20) 
5.416(4) ? 
5.853(15) ? 





TABLE VI. Hyperfine splittings of P levels. 


All values are given in MeV. 




States 


AM D ^AM D AM B 




AAfn ^AMn 

US mh US 


AM Ba 


ip 2 - m 

IPl - lPo 


45 14 14 
63 20 19 




45 14 
61 19 
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FIGURES 



2.6- 



2S 



^ 2S1 

,.-'".'--- 2S0 

-*-'-'"" .... 1P1 

..--'"" . 1P2 

lP(l/2) ---'- ;"-" - IPO 

iM ,.'<:'-—-"" 1P1 

lP(3/2) '-'"" 



2.2- 



1S -c- 

2.0J \ 1S1 



1.8- 



1SO 



mQ— >infinity with 1/mQ corrections 



FIG. 1. The ordering pattern of D meson states. The mass scale is in GeV. 
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2.7- 



2.5- 



2.3- 



1.9- 



.- 2S1 

,----;.'-'-'--- 2S0 

2S "''""" 1P1 

--.- IP? 

lP(l/2) «==:: r . ::;::: - v { v \ 

lP(3/2) '-'-'"'' 1P0 



IS 



2.1-1 \ 1S1 

^ ISO 



mQ— >infinity with 1/mQ corrections 



FIG. 2. The ordering pattern of D s meson states. The mass scale is in GeV. 
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5.8- 



5.6- 



5.2- 



2S1 
2S0 



2S 



lP(l/2) =.-_-_— 1P1 

::: : ipo 

--.-.--= 1P2 

.-===- — 1P1 

lP(3/2) iri 



5.4- is 



1S1 
ISO 



mQ— >infinity with 1/mQ corrections 



FIG. 3. The ordering pattern of B meson states. The mass scale is in GeV. 
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5.9- 



5.7- 



5.3- 



2S1 

__—_-_-----------' 2S0 

2S '"* 

lP(l/2) ^r: : - }g 

...-«=====-"= = IPO 

lP(3/2) — 1P1 



5.5- IS 



1S1 
ISO 



mQ— >infinity with 1/mQ corrections 



FIG. 4. The ordering pattern of B s meson states. The mass scale is in GeV. 



